Astronomy & Astrophysics manuscript no. n4395aaki 


© ESO 2010 


February 10, 2010 





The Suzaku broadband X-ray spectrum of the dwarf Seyfert galaxy 

NGC4395 

K. Iwasawa 1 *, Y. Tanaka 2 , and L.C. Gallo 3 



> 
(N 

\o 

O 

(N 

(N 
O 
O 



INAF-Ossservatorio Astronomico di Bologna, Via Ranzani, 1, 40127 Bologna, Italy 
e-mail: kazushi . iwasawa@oabo . inaf . it 

Max-Planck- Institut fur extraterrestrische Physik, GiessenbachstraBe 1, 85748 Garching, Germany 
e-mail: ytanaka@mpe.mpg . de 

Department of Astronomy and Physics, Saint Mary's University, Halifax, NS B3H 3C3, Canada 
e-mail: lgallo@ap . smu . ca 



Preprint online version: February 10, 2010 



o 



o 

u 

6 
*h 
■i— > 

c3 



ABSTRACT 



We present a Suzaku observation of the dwarf Seyfert galaxy NGC4395 with an estimated black hole mass of ~ 1O 5 M . Rapid and 
strong X-ray variability with an rms amplitude of ~ 60 per cent is observed in the 0.4-10 keV band with the XIS cameras. The shape 
of the light curve appears to depend on energies. The hard X-ray emission is detected up to 35 keV with the HXD-PIN detector 
at a similar flux level as observed with the INTEGRAL IBIS. The X-ray spectrum below 10 keV is strongly absorbed by partially 
ionized ~ 35 erg scnr 1 ) gas with a mean equivalent hydrogen column density of ~ 2 x 10 22 cirT 2 , when a simple absorption model 
is applied. The spectral shape is also strongly variable but not a simple function of the source brightness. The spectral variability 
appears to be accounted for mainly by continuum slope changes, but variability in the ionized absorber may also play some part. The 
apparently flat spectral slope of T =s 1.4 below 10 keV, obtained after correcting for absorption, is marginally inconsistent with the 
r ~ 2 inferred from the 14-35 keV PIN spectrum. If the true spectral slope had been as steep as that measured in the hard X-ray band, 
there would have been an extra absorption component, which we are unable to detect. Combined with the INTEGRAL measurements, 
the hard X-ray emission above 10 keV exceeds the optical emission in terms of luminosity and dominates the broadband energy 
output, unless a large excess of UV disk emission is yet to be detected in the unobservable band. A weak Fe K line is seen at 6.4 keV 
with the average equivalent width of 1 10 eV, which does not show significant flux changes over the 3-day observation. 

Key words. Galaxies: active - Galaxies: individual (NGC 4395) - X-ray: galaxies 



1. Introduction 

NGC4395 is a dwarf galaxy located at the distance of 4.0 Mpc 
(Thim et al 2004) that is known to host a very low luminosity 
Seyfert nucleus (Filippenko, Ho & Sargent 1993). The active nu- 
cleus is located at the centre of the galaxy, which has no obvious 
bulge component. Unlike LINER-like, low-luminosity AGN, the 
optical/UV spectrum of NGC4395 shows broad Balmer emis- 
sion and high-excitation lines, reminiscent of genuine Seyfert 
nuclei but of higher luminosity (e.g., Ho, Filippenko & Sargent 
1997). Using the relation usually applied to higher luminosity 
AGN (Kaspi et al 2000), the bolometric luminosity estimated 
from the 5100 A luminosity is 5.4 x 1 40 erg s _1 (Peterson et al 
2005), which can be powered by a black hole of a few 100 M if 
it accretes close to the Eddington limit. 

There have been indications that the mass of the black hole 
in NGC4395 may be 10 5 M Q or lower, including by X-ray tim- 
ing analysis (Kraemer et al 1999; Lira et al 1999). The mass 
derived from the reverberation mapping technique is found to be 
(3.6+ 1 . 1)X 10 5 M o (Peterson et al 2005). However, the low bulge 
velocity dispersion < 30 km s _1 (Filippenko & Ho 2003) devi- 
ates from the extrapolation of the M-cr relation (e.g., Tremaine 
et al 2002) when this estimate of the black hole mass is adopted. 
Some researchers (e.g., McHardy et al 2006) propose a lower 
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mass, of even an order of magnitude below the value obtained 
from reverberation mapping. 

As pointed out by Peterson et al (2005), the Eddington ratio 
(L bol /L Edd ) is low 1.3 x 10- 3 [M BH /(3.6 x l^Afo)]" 1 . Provided 
that the estimate of the bolometric luminosity derived from the 
5100 A luminosity is correct, even for the lower mass esti- 
mate (~ 1O 4 M ), the accretion rate is the order of 10~ 2 of the 
Eddington limit. In terms of the accretion rate, the X-ray source 
in NGC4395 corresponds to the low/hard state of X-ray binaries. 

The soft X-ray emission observed with ROSAT (Moran et al 
1999, Lira et al 1999) is weak largely because of strong absorp- 
tion (Iwasawa et al 2000). The absorbing material is partially 
ionized (warm absorber) and neutral (cold) absorption is usually 
small (a few times 10 20 ctrT 2 in A^h)- However, an occasional 
increase in the cold absorbing column, perhaps due to a small 
cloud passing across the line of sight, has been observed in one 
of the XMM-Newton observations (Dewangan et al 2009). 

The X-ray emission is known to be extremely variable on 
short timescales. The temporal properties of the X-ray emission 
of NGC4395 was studied in detail using a long XMM-Newton 
observation with a consecutive 90 ks duration (Vaughan et al 
2005). The power spectrum of the flux variation shows a break 
at ~ 2 x 10~ 3 Hz, which is much higher than those for other 
AGN, suggesting a low black hole mass. The rms amplitude of 
the soft X-ray flux variability exceeds 100 per cent. 

In comparison with other AGN, the X-ray spectrum of 
NGC4395 below 10 keV appears to be very hard even after cor- 
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Fig. 1. The X-ray light curve of NGC4395 during the Suzaku observation. The background-subtracted counts in the 0.4-10 keV 
band from the three XISs are summed together. Each timebin is of a 100 s duration. 



recting for absorption, F ~ 1.4 (e.g., Shih et al 2004; Moran et 
al (2005) suggested an even lower value of F ~ 0.6). Despite 
the relatively low 2-10 keV flux (~ 4 x 10" 13 erg s" 1 crcT 2 ), 
NGC4395 was detected by INTEGRAL IBIS (Beckman et al 
2006; Bird et al 2007) in 20-40 keV and 40-100 keV bands. This 
high-energy X-ray detection may imply that the hard spectrum 
continues into the hard band. The measurement of the spectral 
shape above 10 keV, where it is less affected by absorption, is 
important to assess the origin of the hard intrinsic spectrum. We 
observed NGC4395 with Suzaku to measure the spectral shape 
above 10 keV and characterise the flux and spectral variability 
in the softer band by using the X-ray CCD spectrometer over a 
continuous period of more than ~ 2.5 days. 



2. Observation and data reduction 

NGC4395 was observed with Suzaku (Mitsuda et al 2006) on 
2007 June 02 - 05, for a total duration of 230 ks. The XIS cam- 
eras were operated with the Space-Row Charge Injection in the 
standard observation mode. The data were taken by the three 
XISs: two FI CCD cameras (XIS0 and XIS3), one BI CCD cam- 
era (XIS1), and the HXD-PIN. The useful exposures obtained 
from the XIS and the HXD-PIN were 101 ks and 95 ks, re- 
spectively. NGC4395 was observed in the "HXD nominal" posi- 
tion on the detector plane, which is optimised for the maximum 
throughput to the HXD rather than the XIS. 

The data from version 2 processing were used in the anal- 
ysis. The cleaned event files produced by the data processing 
with the standard data selection criteria were used. The non X- 
ray background of the PIN detector was estimated using the 
"tuned" background modeQ The background-corrected count 
rates observed with the three XIS in the 0.4-10 keV band are 
0.107+0.001 (XIS0), 0.116±0.001 (XIS1), 0.109+0.001 (XIS3), 
and 0.031 ± 0.002 with the HXD-PIN in the 14-40 keV band. 



1 http://www.astro.isas.ac.jp/suzaku/process/caveats 

2 http://www.astro.isas.ac.jp/suzaku/analysis/hxd/pinnxb/tuned/ 



3. Results 

3. 1 . X-ray variability in the XIS band 

The X-ray source in NGC4395 is strongly variable during the 
Suzaku observation, as previously observed with other instru- 
ments. Since there are many irregular interruptions in data sam- 
pling because of the low-Earth satellite orbit, the Suzaku data 
are not ideal for a detailed timing study of this rapidly variable 
object on short timescales. We refer to Vaughan et al (2005) 
for such a study with the uninterrupted 90-ks XMM-Newton 
data. Despite the low observation efficiency, the Suzaku data pro- 
vide continuous coverage of the source for more than twice the 
XMM-Newton observation. 

Fig. 1 shows the 0.4-10 keV light curve of NGC4395 ob- 
tained from the three XIS cameras. The resolution is 100 s and 
the count rate is the sum of the three XISs with background cor- 
rection. Large amplitude flux variations are clearly seen. Flux 
doubling in 100-s and shorter intervals is common (Fig. 2) and 
consistent with the characteristic timescale of the variability of a 
few 100 s (Vaughan et al 2005). 

While the fractional variability throughout the observation is 
^var = 0.60 ± 0.06 (see Vaughan et al 2005 for a comparison 
with the XMM-Newton results), the variability over a shorter 
duration can be examined for this rapidly variable object. Fig. 3 
presents the variation in fractional variability amplitude for 5 ks 
intervals in a time sequence obtained from the light curve in Fig. 
1 . A few intervals of excess variability correspond to periods of 
strong, very short-term flaring within 1000-s or shorter duration, 
notably in the hard X-ray band (see Fig. 4). The typical (e.g., me- 
dian) source flux during those periods is low. According to the 
correlation between total rms and flux that is generally found 
in black hole systems (Uttley & McHardy 2001), the quantities 
measured here are on average expected to remain constant. The 
measured values are clearly variable as expected from the ran- 
dom realisation of a source with a red-noise spectrum in flux 
variation (e.g., O'Neill et al 2006). The mean and standard devi- 
ation of F wiI in Fig. 3 is 0.32 and 0.13, respectively, as indicated 
in the figure. Nominal fluctuation of F vax could be represented 
by the standard deviation. All measurements apart from one are 
within 3cr fluctuation. The largest variability recorded in the 80 
ks region (Fig. 2) exceeds 3<x. 
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Fig. 2. The blow-up of part of the light curve in Fig. 1 around 84 
ks where a rapid and large-amplitude flare occurs. The resolution 
of the light curve here is 50 s. This strong flare causes the largest 
F val in Fig. 3. 
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Fig. 3. The variation in the fractional variability amplitude with 
5-ks intervals, obtained from the light curve in Fig. 1. The mean 
value (0.32) and the standard deviation are marked by the dashed 
line and dotted lines, respectively. 
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Fig. 4. The X-ray light curves in four energy bands 0.4-1 keV, 
1-2.5 keV, 2.5-5 keV, and 5-10 keV (from top to bottom) nor- 
malized to the mean count rate. The time resolution is 1 ks. The 
mean count rates in the respective bands are 0.019, 0.109, 0.118, 
and 0.079 count s _1 . The y-axis range for each panel is adjusted 
for visual clarity and varies between the panels. The relative am- 
plitude of the respective light curves can be quantified in Fig. 
5. 



The shape of the light curve depends on the chosen energy 
band. Light curves with 1 ks resolution in the four bands 0.4-1 
keV, 1-2.5 keV, 2.5-5 keV, and 5-10 keV are produced and nor- 
malized to their respective mean count rates (Fig. 4). The dis- 
tribution of the source count rate is far more skewed in the soft 
band than in the hard. This is demonstrated in Fig. 5, in which 
only the 1-2.5 keV and 5-10 keV light curves are displayed for 
clarity. The distribution in the 1-2.5 keV curve is well matched 
to log-normal, while that for the 5-10 keV band is close to being 
a Gaussian. The (rebinned) correlation diagram on the right also 
demonstrates that 1-2.5 keV flux is correlated with 5-10 keV 
flux exponentially rather than linearly, reflecting the difference 
in light-curve shape mentioned above. The 0.4-1 keV band dis- 
tribution is even more skewed than at 1-2.5 keV, but the low 
count rate means that distortion by the Poisson noise could be 
responsible for the extra skewness. 

It is known that the light curves of accreting compact ob- 
jects generally exhibit a log-normal distribution (e.g., Uttley, 
McHardy & Vaughan 2005), which is skewed towards the 
brighter end, as seen in the 1-2.5 keV light curve. We note that 
where a light curve has a log-normal nature, the asymmetry of 
the count distribution will be reduced as the fractional variabil- 
ity becomes significantly small (see simulated light curves with 
the fractional variability of 17, 41, and 64 per cent in Fig. 4 in 
Uttley et al 2005). However, in NGC4395, the 5-10 keV band 



variability remains strong, and the reduction in the fractional 
variability (from 70 per cent in the soft band to 40 per cent in 
the 5-10 keV, Fig. 6) may not be large enough to account for 
the change observed. As described below, variability in the ion- 
ized absorber could alter the shape of the distribution. Another 
cautionary note is whether the observation is sufficiently long to 
cover the timescale on which the stationarity of the time series 
is achieved, unlike measurements for Galactic binaries. 

The log-normal shape of a distribution is expected in multi- 
plicative products of random events, and, in the context of black 
hole systems, a model in which fluctuations in the accretion rate 
are amplified as they propagate inwards has been proposed (e.g., 
Lyubarskii 1997). In contrast, the symmetric shape of the hard- 
band light curve may be unusual. A superposition of random 
shots is expected to produce a Gaussian distribution, as seen in 
the hard X-ray light curve. While the soft and hard X-ray emis- 
sion in NGC4395 appear to behave independently on some oc- 
casions (Fig. 4), they are statistically correlated (Fig. 5 right). 
The correlation occurs in a way that the soft X-ray variability is 
stretched to higher fluxes, i.e., in an exponential manner relative 
to the hard X-ray (Fig. 5 right). Seyfert galaxies generally shows 
an approximately linear correlation in two bands with an offset, 
which is attributed to a constant component such as reflection, 
while a power-law relation, which is generally expected (Uttley 
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Fig. 5. Left: The distribution of the normalized count rates in the 
1-2.5 keV (dash-dotted line) and 5-10 keV (solid line) bands in 
Fig. 3. It is clear that the distribution in the soft-band light curve 
is strongly skewed and well matched by a log-normal shape, 
while that for the hard band is close to being symmetric (i.e., 
a Gaussian). Right: The correlation diagram between the count 
rate in the 5-10 keV and 1-2.5 keV band. The data have been in- 
crementally rebinned according to the 5-10 keV count rate, and 
the five bins contain 27, 27, 27, 27, and 28 measurements in as- 
cending order. Note that the y-axis (1-2.5 keV count rate) is in 
a logarithmic scale while the x-axis (5-10 keV count rate) is in 
linear scale. 
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Fig. 6. The fractional rms variability amplitude (F var ) as a func- 
tion of energy, obtained from the light curves in 1 ks resolution 
(see Fig. 3). 



et al 2005), is more successful in reproducing data for the highly 
variable NGC405 1 (Taylor, Uttley & McHardy 2003). 

We note that in Fig. 4, the two relatively broad, intervals con- 
taining high flux at 50-90 ks and 160-190 ks have different en- 
ergy dependences (also see Section 3.4). Very short-period/rapid 
flickering at the 80-90 ks intervals (where the largest variability 
is recorded in Fig. 3; see also Fig. 2) appears more pronounced 
at higher energies. The flaring at 90-100 ks is only clear for the 
two harder bands. 

The overall variability amplitude is larger at lower energies, 
as can be seen in Fig. 6. The fractional variability is similar 
within the two low energy bands at ~ 70 per cent, and declines 
towards higher energies to ~ 40 per cent in the 5-10 keV band. 

3.2. PIN band variability 

The entire observation is divided into 5 time-segments: the first 
four segments are of 50 ks duration and the last of 30 ks. The 
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Fig. 7. The light curves of NGC4395 in the 13-33 keV band from 
the HXD-PIN (upper panel) and the 9-11 keV band from the 
XIS0 and the XIS3 (bottom panel). The dotted line indicates the 
mean counts over the same interval as that of the HXD-PIN light 
curve. 



13-33 keV count rates from the PIN for these five intervals are 
plotted in Fig 7. The error bars in the count rate are statistical 
only. There appears to be a flux increase by ~ 0.02 ct s _1 in the 
final 80 ks of the observation. In this energy band, 5 per cent 
of the non X-ray background of the PIN is about 0.018 ct s -1 , 
which is the conservative estimate of the uncertainty given the 
exposure of the individual data points. The possible flux vari- 
ation is comparable to this, but the reality of the apparent flux 
variation is uncertain. As a comparison, the light curve for the 
9-11 keV band, which was obtained by the XIS0 and XIS3 (the 
signal from the FI-CCD of the XIS 1 is very noisy in this band), 
is also shown in Fig. 7. Flux variability of 30 ± 5 per cent is 
still present at the highest energy end of the XIS bandpass when 
viewed on the 5 ks time resolution. On the 50 ks timescale, the 
amplitude is lower, but a similar variability pattern between the 
two light curves indicates that the hard X-ray variability in the 
PIN band may be real. 

3.3. The time-averaged energy spectrum 
3.3.1. The XIS data 

A simple power law is unacceptable in describing the full-band 
data mainly because of the strong absorption at low energies 
(Fig. 8). The best-fit model photon index F ~ 0.7 indicates that 
the overall spectrum is unusually hard. Because of the strong 
absorption and hard spectrum, the spectral analysis here focuses 
mainly on the two FI CCD cameras, XIS0 and XIS3. 

Introducing a warm absorber and a Fe K emission-line pro- 
vides a good description of the 0.4-10 keV band spectrum 
(Fig. 9). The power-law continuum has a photon-index of T = 
1.40+°;^. One of Tim Kallman's XSTAR grid tables (grid25, 
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with a Gaussian profile, has a centroid energy at 6.36 + 0.04 keV, 
a line width <x = 0.09 + 0.05 keV (or FWHM- 10, 000 km s 1 ), 
and line intensity 6.4^ ? x 10~ 6 phs 1 cm -2 , which corresponds 
to an equivalent width of EW = 115^2 eV. The quality of fit 
(X 2 = 406.0 for 407 degrees of freedom) is good (Table 1). We 
note that there is a hint of systematic hardening above 9 keV 
(see Fig. 8 and the residual plot in Fig. 9). However, the data in 
this range have small statistical weight in the full-band spectrum 
and the deviation does not have an impact on the quality of the 
overall fit. 

The observed flux is 5.1 X 10~ 13 erg s _1 cirT 2 in the 0.5-2 keV 
band and 4.4 x 10~ 12 erg s" 1 cm" 2 in the 2-10 keV band. The 2- 
10 keV flux is similar to that during the long ASCA observation 
(4.7 x 10~ 12 erg s _1 cirT 2 , Shih et al 2002) and slightly (~ 20 per 
cent) lower than that during the long XMM-Newton observation 
(5.6 X 10~ 12 ergs _1 cm" 2 , Vaughan et al 2005). 



Fig. 8. The 0.5-10 keV spectrum of NGC4395 observed with 
the Suzaku XIS. The data are obtained from the XIS0 and XIS3, 
and plotted in units of flux density after correcting for the detec- 
tor response. This correction may have some uncertainty but is 
accurate enough to show the major spectral features in the spec- 
trum. 
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Fig. 9. The 0.5-35 keV spectral data of NGC4395 obtained 
with the XIS and HXD-PIN. The solid-line histogram shows the 
model of a power law with the Fe K line modified by the warm 
absorber folded through the detector responses, which fits the 
XIS data. The cross-calibration constant is assumed to be 1:1.18 
(XIS:PIN). 



which is the photoionization grid computed by XSTAR with a 
power-law ionizing source and the turbulent velocity of the ion- 
ized gas of 200 km s'nis used to model the warm absorber. 

' . . . ... . ... _ . . <» 1 



The absorbing column of the warm absorber is /Vh,w = 2.1 



+0.3 
-0.2 



10 cm , and the ionization parameter log £ = 1.53+qj (erg cm 

s _1 ). A small cold absorption of /Vh= 5.4x 10 20 cirT 2 , in addition 
to the Galactic column (N H = 1.9 x 10 20 cnT 2 , for the LAB map 
by Kalberla et al 2005) is also required. The Fe K line, modelled 



3.3.2. PIN data 

The 15-35 keV spectrum measured with the HXD-PIN appears 
to be much steeper than that below 10 keV, albeit with a large 
uncertainty. Fitting a power law provides a reasonable fit to the 
data with a photon index of T = 2.2 ± 0.5 Oc 2 = 14.5 for the 8 
degrees of freedom). The inclusion of any spectral curvature or 
break does not improve the fit. The 20-40 keV flux, estimated by 
the simple power-law model, is 5.9 x 10 _12 ergs _1 cirT 2 , which 
is similar to those measured previously by the INTEGRAL IBIS 
(Beckman et al 2006; Bird et al 2007). 

The broad-band spectrum is examined by combining the XIS 
and the PIN spectrum. Since the observation efficiency of the 
PIN is slightly lower than that of the XIS, good time intervals 
for the PIN are also used to produce the XIS spectrum analysed 
here (~ 85 per cent of the original XIS exposure). 

When extrapolating the best-fit model of the XIS spectrum 
(Fig. 9), the lower end of the PIN data are well above the extrap- 
olation. The cross-calibration constant of l.lS^between the PIN 
and the XIS is adopted, which was measured for the Crab neb- 
ula at the HXD nominal pointing position. Provided this cross 
calibration is reliable, a spectral hump in the 10-25 keV range is 
required, which may be supported by the possible spectral hard- 
ening seen at the highest energy end of the XIS data (Figs. 8 and 
9). However, we caution that the uncertainty in the PIN back- 
ground estimate may significantly influence whether we detect 
a spectral hump for this faint source. The source fraction of the 
total PIN counts in the energy band of interest is only ~ 9 per 
cent. +3 per cent uncertainty in the background estimate trans- 
lates into the cross-calibration constant in the range of 0.8-1.6 
being allowed. A close match between the extrapolated model 
and the 15-20 keV PIN data can indeed be obtained when the 
cross-calibration constant is increased to ~ 1.5. 



3.4. Spectral variations 

As inferred by the study of light curves in various energy bands, 
the spectral shape of the X-ray source in NGC4395 is also vari- 
able. Furthermore, the spectral variability does not appear to cor- 
relate with the source brightness (see Section 3.1). 

The throughput of the Suzaku XIS is insufficient to track the 
spectral changes on the time-scale of the most rapid flux vari- 
ability (of the order of 100 s). Instead, the spectral variability 



http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html 



4 ftp://legacy.gsfc.nasa.gov/suzaku/doc/xrt/suzakumemo-2008- 
06.pdf 
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Table 1. Spectral fit to the Suzaku XIS time-averaged spectrum. 

Data T N H iV H ,w f ^FeK cr / FcK AT 2 /dof 

crrr 2 cm -2 ergs~'cm keV keV phs~'crrr 2 

Total 1.40^ 5.4! 2 *xl0 20 2.1^ xlO 22 34^ 6.36^ 0.09^ 6.4^ x 10" 6 406.0/407 

Note: The spectrum was obtained by combining the XIS0 and XIS3 data. The model consists of a power law with a Gaussian for the Fe Ka 
line, modified by warm and cold absorption. The cold absorption column is denoted by Nu, which is in excess of the Galactic column of 
1.85 x 10 20 cm~ 2 . Tim Kallman's XSTAR grid tables ("grid 25") with the turbulent velocity of 200 km s -1 is adopted for the warm absorber 
modelling. 



Table 2. The definition of the four X-ray colours. 
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keV 


keV 




HR1 


0.4-1 


1-2.5 


42.2/45 


HR2 


1-2.5 


2.5-5 


323/45 


HR3 


2.5-5 


5-10 


115/45 


HR4 


0.4-1 


5-10 


117/45 



Note: Each colour indicator (HR) is defined by (H-S)/(H+S), where 
H and S are count rates measured in the two energy bands given 
above, x 2 values for a constant hypothesis are also included. 



was investigated on a timescale of 5-ks using X-ray colours. The 
four energy bands 0.4-1 keV, 1-2.5 keV, 2.5-5 keV, and 5-10 
keV were used in the X-ray colour analysis. The X-ray colour 
is defined to be HR = (H-S)/(H+S), where H and S are the 
background-corrected count rates collected from the three XIS 
cameras in two bands of choice and H is always of the higher 
energy. Four HR (HR1, HR2, HR3, and HR4) are derived as 
shown in Table 2. 

The X-ray colour measurements are plotted in Fig. 10, along 
with the full-band light curve. Four correlation-diagrams be- 
tween the X-ray colours and the source brightness are presented 
in Fig. 1 1 . Overall, there is no clear correlation between the spec- 
tral variability and source flux, although some systematic trends 
can be found: 

1) HR values increase towards lower flux down to the 0.4- 
10 keV count rate of ~ 0.2 ct s" 1 , except for HR1, 2) Below 
this flux, HR drops, except for HR3. The former trend is con- 
sistent with spectral pivoting, i.e., the spectral slope change as a 
function of source brightness. However, the turnover of the HR 
value in the low flux range suggests further complication to a 
simple spectral pivoting. Below 1 keV, there is an emission com- 
ponent from photoionized gas, detected through XMM-Newton 
RGS measurement of OVII (Bianchi & Guaunazzi 2007), which 
is detached from the variable nuclear emission and expected to 
be stable. This could partly explain the HR turnover, namely in 
HR1 and HR4, as the stable component becomes visible when 
the nuclear flux is low. 

In the time sequence, there is considerable flaring in the 40- 
80 ks interval, where the spectrum remains soft, as shown by 
generally low values of HR (Fig. 10). However, the softest pe- 
riod, indicated by the HR curves, seems to last for a while (~ 20 
ks) after the flare peak. Similar HR behaviour may also occur 
during the decay of the large flare at the end of the observation. 
These are seen in HR2, which shows the most significant varia- 
tion (see Table 2), and possibly in HR4. The colour HR2 is taken 
from the energy range where the source count rate is the highest, 
which is probably the main reason for the large significance, but 



Table 3. The XIS 1-10 keV spectra of the six time-intervals. 



Interval 


T 


log? 




r 2-10 


1 


1.50 ±0.07 


1.6 ±0.1 


1.6 ±0.2 


4.8 


2 


1.67 ±0.06 


1.7 ±0.1 


1.2 ±0.2 


5.2 


3 


1.32 ±0.08 


1.3 ±0.1 


2.0 ±0.2 


3.3 


4 


1.08 ±0.06 


1.0 ±0.1 


2.4 ±0.2 


2.6 


5 


1.16 ±0.06 


1.0 ±0.1 


2.4 ±0.2 


4.9 


6 


1.58 ±0.06 


1.6 ±0.1 


1.3 ±0.2 


5.6 


r 




1.47 ±0.04 


1.40 ±0.04 




log^ 


1.5 ±0.1 




1.2 ±0.1 




K 


2.2 ±0.3 


2.2 ±0.1 






X 2 /dof 


671.8/676 


707.2/676 


698.1/676 





Note: Results of fitting a single warm absorber model to the 1-10 
keV XIS spectra from the six time-intervals are presented. One of 
the three parameters, T, and JVw> is a f fee parameter when fitting 
the six spectra, besides the normalizations of the power-law con- 
tinuum. The other two are tied together between the spectra and 
reported below. Observed 2-10 keV flux is also given. 

a In unit of erg s -1 cm. 

b In unit of 10 22 cm" 2 . 

b The 2-10 keV flux in unit of 10~ 12 erg cur 2 s _1 . 



Table 4. The XIS spectra fromthe second and fifth intervals. 



Interval 


r 


logf 






2 


1.75 ±0.07 


1.6 ±0.5 


1.1 ±0.2 




5 


1.23 ±0.07 


1.0 ±0.1 


2.4 ±0.3 




r 




1.56 ±0.09 


1.47 ±0.09 




io g r 


1.4 ±0.3 




1.1 ±0.2 




7V b 


2.2 ±0.5 


2.3 ±0.3 






* 2 /dof 


281.8/309 


293.2/309 


283.4/309 




Interval 


r 


logf 




* 2 /dof 


2 


1.65 ±0.08 


1.5 ±0.2 


2.2 ±0.4 


275.0/308 


5 


1.34 ±0.10 


1.0 ±0.2 


as above 




2 


1.66 ±0.09 


1.3 ±0.2 


1.8 ±0.2 


275.5/308 


5 


1.34 ±0.09 


as above 


2.3 ±0.3 





Note: Results of joint model fitting to the XIS spectra from the sec- 
ond and fifth intervals are shown. For the spectral fitting, the proce- 
dure is the same as that described in Table 3. The lower part shows 
the results when the spectral slope and the warm absorber parame- 
ters (Nw, log £) are left to vary, but either of the two warm absorber 
parameters is tied between the two spectra. 

a In unit of erg s -1 cm. 

b In unit of 10 22 cm" 2 . 
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Fig. 10. The temporal variations in the 0.4-10 keV count rate 
and the four X-ray colours (HR1, HR2, HR3, and HR4) at 5 ks 
resolution (from top to bottom). 



where also the warm absorber modifies the spectrum the most. 
In the simplest picture of photoionization, the ionization param- 
eter should change in response to variations in the illuminating 
source and the ionic column of the absorber should also change 
on the recombination timescale. The complicated HR behaviour 
may thus also originate from the variability in absorption. 

In an attempt to identify the main driver of the spectral vari- 
ability in NGC4395, we divided the observation into six inter- 
vals and examined their spectra, by applying the simple model 
which was used to model the time-averaged spectrum (Table 1). 
All the time intervals are of a 40 ks duration each, except for the 
last interval, which has a 30 ks duration. The second and fifth 
intervals roughly correspond to the periods, where the spectrum 
is persistently soft and hard, respectively, according to HR2. 

The six spectra were then fitted jointly by restricting the 
choice of the model free parameters to the power-law normal- 
isations and one parameter of choice among photon index T, 
the ionization parameter and the absorption column density 
Nw of the warm absorber, either of which could be the driver of 
the spectral variability. To avoid possible contamination by pho- 
toionized gas at the low flux intervals, the data below 1 keV were 
excluded. The cold absorption was fixed to the best-fit model 
value for the time-averaged spectrum (see Table 1), and all other 
parameters are tied together between the six spectra. Results are 
summarised in Table 3. The highest quality of fit is obtained 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 
0.4-1 OkeV (ct/s) 0.4-1 OkeV (ct/s) 

Fig. 11. The three X-ray colour indicators HR1, HR2, HR3, and 
HR4 (see Table 2) plotted against the 0.4-10 keV count rate. 



when T is a free parameter, suggesting that the slope change 
might be the primary driver of the spectral variability. 

However, as discussed for the HR analysis, the slope change 
may not be the unique cause of the spectral variability. We then 
took the second and fifth intervals, which have almost identical 
5-10 keV flux, =! 2.9 x 10~ 12 erg s -1 cirT 2 (see Fig 4) but signif- 
icantly different values in HR2 (see Fig. 10). The same exercise 
performed for these two intervals infers that a change in the con- 
tinuum slope and the absorption depth are comparably good so- 
lutions (see Table 4). When either £ or Af w were allowed to vary 
between a pair of intervals, in addition to T, the quality of the 
fit improved by a similar degree. The reduction in x 2 relative to 
the fit, where only F is variable between the two spectra (Table 
4), is Ax 2 ~ -6, for which a F-test infers 99 per cent confidence 
in the one additional term. No additional improvement in the fit 
was attainable when both ^ and Nw were allowed to vary. This 
may indicate that changes in the warm absorber, at least in one 
of the two parameters, could be a secondary source of the spec- 
tral variability. Which parameter varies of either £ or A^w in the 
warm absorber is difficult to determine, because their features 
are degraded at the CCD resolution. The higher value of £ for 
the spectrum with higher flux is qualitatively consistent with the 
photoionization of the absorbing medium. 

Finally, we comment on the spectrum at the minimum of 
the light curve at ~ 150 ks. During this flux minimum, HR1 
shows spectral softening, while HR3 shows strong hardening. 
The softening could be explained by the increased relative con- 
tribution of the stable, photoionized-gas emission to the soft 
band. On the other hand, the larger value of HR3 is caused by 
rapid flaring, which only occurs in the 5-10 keV band (see Fig. 
4). Naturally, this hard X-ray flare has an extremely flat spec- 
trum (r = 0.25 ± 0.15), which makes it difficult to explain using 
the absorption model. In this brief period, the variable contin- 
uum does not switch off completely, as significant variability is 
still present (Fig. 2). Therefore, explanation of the hardening in 
HR3, invoking a constant component (in this case, the cold re- 
flection from distant matter), which would appear mainly in the 
5-10 keV band, does not apply. 
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Fig. 12. The count rate spectra for the two intervals (40-80 ks 
and 160-200 ks), where the 5-10 keV flux is comparable but the 
HR2 shows largely different values. These two intervals corre- 
spond to the second and fifth intervals in Table 3, and investi- 
gated in Table 4. The spectra are obtained by combining the data 
from the XIS0 and XIS3. The solid-line histograms indicate the 
best-fit model in which the power-law normalization and the ab- 
sorption column (AT W , see Table 4) are the only independent pa- 
rameters between the two. The bottom panel shows the residuals 
in the form of the ratio of the data to the model for the two spec- 
tra. 



4. Iron line variability 

The iron K line variability was investigated in time and also at 
different continuum flux levels. Since the line is weak (Table 
1), it is not possible to track the response of the line flux to the 
extremely rapid continuum variation. The XIS spectra from three 
time intervals i-1 of 0-90 ks, i-2 of 90-160 ks, and i-3 of 160-230 
ks were extracted and the line flux measured. The intervals, i-1 
and i-3 contain periods of strong flaring, while i-2 represents 
relatively quiet periods. 

The spectral data are restricted to the 3-10 keV range and the 
continuum and the iron line are modelled by an absorbed power 
law and a Gaussian, respectively. The line flux is not corrected 
for absorption. The line fluxes in the three intervals agree and 
show no sign of variability, as shown in Table 5. 

Next, the source light curve was sliced into three count-rate 
ranges and the line flux from each source brightness range was 
measured. The count-rate slicing was applied to the light curve 
in Fig. 1 (in the 0.4-10 keV band with the 100-s time resolu- 
tion) with the count rate ranges ofL: < 0.30 cts _1 ,M: 0.30-0.55 
cts -1 , and H: > 0.55 cts -1 . The line flux measurements were 
performed in the same way as for the time-ordered spectra. The 
results are presented in Table 5. There is a hint of line flux in- 
crease in the data at the highest source brightness, but within 
the uncertainty, the line flux does not show significant variabil- 
ity over factor of ~ 3 change in average flux. At least, a linear 
correlation between the continuum and line flux can be ruled out. 
Table 3 also shows exposure time, the collected counts, and the 
flux in the 3-10 keV band. 



o 




-I - 

O r 



10 9 10 10 10" 10 12 10 13 10 14 10 15 10 16 10 17 10 18 10 19 
Frequency (Hz) 

Fig. 13. The energy distribution of NGC4395 from radio to 
hard X-rays. The X-ray data (open squares) are measured with 
Suzaku. The dotted histograms indicate the spectra during the 
low- and high-flux states observed in the Suzaku observation. 
The 20-40 keV and 40-100 keV data points (open circles) are 
from INTEGRAL IBIS (Bird et al 2007). The other points are 
taken from NED, selecting photometric data from the smallest 
aperture: SDSS (u, g, r, i, z), HST/ACS (F330W filter, 0.3 arcsec 
aperture, Munoz Martin et al 2007), 2MASS (J, H, K s , 14 arcsec 
aperture), MMT (N, 5.3 arcsec aperture, Maiolino et al 1995), 
VLBA (5 GHz, Nagar, Falcke & Wilson 2005), and FIRST (1.4 
GHz, Becker, White & Helfand 1995). IRAS fluxes are not in- 
cluded since they are believed to originate in the galaxy, not from 
the nucleus. The energy distribution of NGC4395 is unusual for 
Seyfert galaxies or quasars (e.g., Elvis et al 1994), the energy 
output at 10 keV or above being comparable to that in the opti- 
cal band. 

Ideally, the flux slicing should be performed using the ioniz- 
ing flux, e.g., the 7-11 keV band flux responsible for ionizing Fe. 
The much lower count rate in this hard band means however that 
the time resolution becomes poorer. We made three flux slices 
based on the 7-11 keV light curve with 10-ks resolution, and in- 
vestigated the Fe line flux. While the dynamic range of the three 
flux slices are only a factor of 2, the line flux remains constant 
between the flux slices, confirming the result obtained above. 
The iron K line is resolved and the inferred FWHM (~ 10000 km 
s _1 ) implies a smaller line-emitting radius than the BLR (FWHI 
5000 km s" 1 , Ho et al 1997). The lack of significant variability 
in the line flux is unexpected in this compact system. Additional 
studies of a wider range of ionizing flux variations will be needed 
to verify the result. 

5. Discussion 

5.1. Origin of spectral variability 

Strong spectral variability is observed in the XIS band. The gen- 
eral trend is that the spectrum is steeper when the source is 
brighter, and vice versa, although a close inspection shows that 
it is not a simple function of the source brightness (Fig. 11 and 
Table 3). The spectral analysis of the six time intervals (Table 
3) appears to show that the spectral variability can be largely ac- 
counted for by continuum slope changes. Variability in the warm 
absorber is less likely, but a good alternative source of the vari- 
ability and may also change the spectrum. However, for the rapid 



K. Iwasawa et al.: The Suzaku broadband X-ray spectrum of the dwarf Seyfert galaxy NGC4395 



9 



Table 5. Iron K line variability. 



Data Range 


Exposure 


Counts 




(1) 


(2) 


(3) 


(4) 


TIME ORDER 








i-1 0-90 


40.0 


5243 


4.0 


i-2 90-160 


29.8 


2450 


2.5 


i-3 160-230 


29.0 


4117 


4.3 


FLUX ORDER 








L < 0.30 


51.8 


3592 


2.2 


M 0.30-0.55 


31.7 


4772 


4.4 


H > 0.55 


16.2 


3681 


6.7 



^FcK 

(5) 

7 5+ 2 - 3 
8.l5| 



7 + ' : ^ 
-M 
9.4+** 



(1) in unit of 10 3 s in the light curve in Fig. 1 for the time intervals, 
or ct s~'for the source count rate range 

(2) 10 3 s 

(3) background-corrected integrated counts 

(4) 3-10 keV flux in 10^ 12 ° 



erg s 



cm • 

-2 



(5) Fe K line flux in 10 6 ph s ' cm 
Note: The error in the line flux is of the order of Icr. 



changes observed in NGC4395, it is generally difficult to draw 
strong conclusions using spectra of the CCD spectrometer's res- 
olution as to whether the warm absorber is actually responsible 
for the variability. 

The continuum slope can vary rapidly, as it is probably con- 
trolled by the heating/cooling in the hot corona at small radii 
of the accretion disk. Any change in the physical conditions 
of the corona can propagate approximately in a light crossing 
time, which is a few tens of seconds for the radius of 10 grav- 
itational radius in NGC 4395. On the other hand, the warm ab- 
sorber would need a longer time, the order of the recombination 
timescale, which depends on the density of the absorbing mat- 
ter, to respond to the variability of the illuminating source. The 
X-ray colour HR2, derived from the 1-2.5 keV and 2.5-5 keV 
data, would be sensitive to any changes in the warm absorber. 
Possible delays in the response of the hardness ratio to the source 
flux variation seen in some occasions may be a signature of vari- 
ations in the warm absrober. 

we now consider that this possible delay of 2 x 10 4 s is 
the response time of the warm absorber. For an ionizing lumi- 
nosity of 3 x 10 4() ergs -1 , the calculation of the recombination 
timescales of the relevant photoionized gas, using the formulae 
given in Otani et al (1996) and Krolik & Kriss (2001), yields the 
density of n ~ 10 7 cm" 3 and the distance of the absorbing matter 
r ~ 3 x 10~ 3 pc. This density is relatively high but lower than that 
inferred for the variable warm absorber in NGC 4051 (Krongold 
et al 2007) in which the X-ray spectral variability is attributed to 
variable absorption (but see Ponti et al 2006). 

5.2. Hard X-ray spectrum and broad-band energy distribution 

The hard X-ray spectrum of NGC4395 above 10 keV, measured 
with the HXD-PIN, is F ~ 2, the typical value for Seyfert 
galaxies, while the spectrum below 10 keV is much harder; 
r =! 1.4, even after correcting for the absorption (Table 1). 
Although the PIN slope is not well constrained, this might be 
indicative of a spectral break around 10-20 keV or an extra ab- 
sorption component that we failed to model. A comparison with 
two INTEGRAL data points in the 20-40 keV and 40-100 keV 
bands (Fig. 13) suggests that the hard X-ray spectrum may actu- 
ally continue towards 100 keV with a slope r ~ 1.5, similar to 
that obtained from the XIS data. However, since the INTEGRAL 



data were not taken contemporaneously, this is still uncertain due 
to possible flux variability. A stronger constraint on the hard X- 
ray spectral shape would therefore be desirable. NGC 4395 is 
not in the Fermi LAT Bright Source List |^] which at least rules 
out such a hard spectrum not extending to 100 MeV. 

The multiwavelength spectral energy distribution (SED) is 
constructed by taking photometric data collected from NED 
(Fig. 13; see also Moran et al 1999; Lira et al 1999; Iwasawa 
et al 2000). Only the data taken from a small aperture were used 
to minimize the contamination from the host galaxy, which can 
be large for this low luminosity nucleus. The IRAS data points 
are discarded since most of the IRAS fluxes are believed to orig- 
inate in the galaxy (Lira et al 1999). 

Moran et al (1999) pointed out the anomalous SED shape 
of NGC4395, which differs from that for the luminous QSOs 
(e.g., Elvis et al 1994) and also from the low luminosity AGN 
(Ho 1999), as suggested by the absence of an optical-UV bump 
and the radio quiet nature of the galaxy (the radio luminosity 5 
orders of magnitude below the optical; see also Wrobel & Ho 
2006). Dust reddening, which might suppress the optical light, 
was discussed by Moran et al (1999), who did not find com- 
pelling evidence of significant reddening. A remarkable feature 
of the SED in Fig. 13 is that the peak of the energy output oc- 
curs at the hard X-ray range, unlike normal Seyfert galaxies and 
QSOs, in which optical-UV emission dominates the bolomet- 
ric output. Without soft seed photons from the accretion disk, it 
would be difficult to explain the production of this strong hard 
X-ray emission with the thermal Comptonization model (e.g., 
Haardt & Maraschi 1991), unless the disk thermal emission es- 
capes our detection into the unobservable shorter wavelength 
UV range, given the small black hole mass, e.g., the blackbody 
emission peak would move by a half decade in wavelength, rel- 
ative to Seyfert galaxies with a black hole mass of ~ 1O 7 M . 

The distinctive SED shape of NGC4395 implies that the use 
of the empirical relation to estimate the bolometric luminosity 
from the 5100 A luminosity, which is calibrated by the data 
for luminous AGN, may be inappropriate. The estimated 0.5- 
100 keV luminosity is ~ 4 x 10 40 erg s _1 , which happens to be 
comparable to the empirical estimate of the bolometric luminos- 
ity. If the bulk of the bolometric luminosity indeed comes from 
the hard X-ray emission, the Eddington ratio remains low, 10~ 3 - 
10~ 2 , depending on the Mbh estimate, and lies in an interme- 
diate range between the low accretion-rate (the Eddington ratio 
of < 10~ 3 ) objects such as M87 and normal Seyfert/QSOs. The 
radio quietness of the nucleus suggests that the accretion flow 
is probably not in the form of an advection-dominated type as 
proposed for objects with a much lower accretion rate. However, 
we note that NGC4395 is located on the "fundamental plane of 
black hole activity" proposed by Merloni, Heinz & Di Matteo 
(2003). The fundamental plane is parameterised by black hole 
mass, radio luminosity, and X-ray luminosity. Fitting this rela- 
tionship suggests that the dim optical emission could be an un- 
usual feature for this active nucleus. 
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